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Abstract 



A Schwinger-Dyson equation for the quark propagator is derived in the 
context of a Bethe-Salpeter second order formaUsm developped in preceding 
papers and of the Minimal Area Law model for the evaluation of the Wilson 
loop. We discuss how the equal time straight line approximation has to be 
modified to include correctly trajectories going backwards in time. We also 
show, by an appropriate selection of the solution of the SD equation, that in 
the limit of zero quark mass chiral symmetry breaking and a zero mass pseu- 
doscalar meson actually occur. The inclusion of backward quark trajectories 
proves to be essential to make the model consistent with Goldstone's theorem. 
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I. INTRODUCTION 



The problem of the chiral symmetry breaking in QCD [|I| and that of the consistency 
with Goldstone's theorem of approximations in the kernels of the Schwinger-Dyson (SD) and 
Bethe-Salpeter (BS) equations has been discussed by various authors |@J^. To my knoledge, 
however, in all such papers only the perturbative part of the kernel was actually derived from 
QCD. To this an ad hoc infrared singular term was added only for the sake of convenience 
and to make the treatment of the resulting equations as easy as possible. 

On the other hand, taking advantage of a four-dimensional path integral representation 
for a second orc/er quark-antiquark Green function H{xi, x^] Hi, 1/2), in a preceding paper we 
succeeded in obtaining a BS like equation for such quantity entirely from first principles, 
a part the use of the so called Minimal Area Law (MAL) model for the evaluation of the 
Wilson loop |J|. 

In this paper I first show that along similar lines a Schwinger-Dyson equation can be 
derived for the colorless second order quark propagator H{x — y) occurring in the BS equa- 
tion. Then I discuss how the equal time straight line approximation to the minimal area 
has to be modified in order to make the formalism consistent with Goldstone's theorem in 
the limit of zero quark masses. 

The mentioned path integral representation for H{xi,X2', 2/2) is in turn a consequence 
of the similar Feynmann-Schwinger (FS) representation for a one particle propagator in an 
external field. Its interest lies in the fact that the gauge field appears in it only in terms of 
the Wilson loop made by the quark and the antiquark trajectories (connecting yi to xi and 
y2 to X2) closed by two Schwinger strings 

W = -(TrPexp{^£/ i" dx^A^}) . (1.1) 



As usual in ( |1 . 1| ) = |v4j^A° is a colour matrix, P the ordering prescription along the loop 



and the expectation value stands for a functional integration on the gauge field alone. The 
MAL model consists in assuming that the logaritm of the Wilson loop can be written as the 
sum of a perturbative and an area term 

ilnW = t{\nW)pen + aSrain, (1.2) 

a being the so called string tension and S'min the minimal area enclosed by Tgg. In practice 
5*111111 is replaced by its equal time straight line approximation, consisting in setting 

'S'min = dTi dT25{zio ~ 2;2o)|zi - Z2I dX^^zlf^z^o - (Ai'iTi2o + (1 - A)i'2T2;io) , 

(1.3) 

where Zi = Zi{ti) and Z2 = Z2{t2) (with < Ti < si and < r2 < S2) are the parametric 
equations for the quark and the antiquark trajectories respectively and Zit, i*2T stand for 
the transverse components of Zi and Z2 [i^x = [5^^ — f^f-^^Zj with r = Zi — 2*2 and = 

In the limits X2 —>■ Xi and y2 yi (i.e. when the initial and respectively the final points 
of the quark and the antiquark trajectories coincide), Eq. ( |1.3| ) is correct only up to the 
second order terms in Zi and Z2 in the general case, but it becomes exact in two different 
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significant situations: when the two trajectories he on a plane (with one time and one space 
dimension) and when they form a double elicoid. This fact may justify its use even in a 
fully relativistic case. 

Notice however that the right hand side of ( |1.3|) is of a purely qq interaction type, while 
the perturbative contribution contains both an qq interaction term and two corresponding 
q and q selfinteraction terms. As a consequence of this lack of symmetry, in the limit of 
zero quark mass ( [L.3D proves to be inconsistent with the Goldstone theorem, if dynamical 
symmetry breaking occurs. We shall see that we obtain consistency if we replace (|1.3|) with 

Smin = dn dT26{zio - 22o)|^i - ^2 1 e(i;io)e(i2o) 
Jo Jo 



X ^ dX^zf^z^o - iXznZ2o + (1 - X)z2TZiofY " 

/ dn dT[6{zio - z[q)\zi - z'^\e{zio)e{z[Q) 
Jo Jo 

X f (iAjiioifo - (Afixiio + (1 - >')ziTZioy}^ - 



S2 rT2 

dT2 / dT^5{z2o - 4o)l^2 - 4|e(2^2o)e(4o) 



JO 



X 



J dX{z2QZ2Q - {Xz2Tz'2o + (1 - A) ^2t2^2o) ^ } ' , 

(1.4) 

where we have used z[ and 4 for zi{t[) and Z2{t2) respectively and e{t) for the sign function. 



In fact ( |1.3| ) breaks down necessarily if the trajectories can go backwards in time. On 
the contrary, if they never do that, (|1.4| ) becomes identical to (|1.3| ) (since in this case the 
first term is positive and the other two vanish), but for a plane or an elicoid loop it remains 
exact, as we shall see, even if they do. The correct inclusion of such more general trajectories 
turns out to be essential and affects substantially the SD equation. 

The resulting SD equation turns out to be equivalent to a system of four non linear inte- 
gral equations envolving four scalar quantities, hQ{k), . . . h^i^k), which appear as coefficients 
of product of Dirac matrices with different tensorial character in the expression of the quark 
propagator. We are not able to handle rigorously such a system. However, its struture 
suggests the possibility of four different classes of solutions corresponding to different sets 
of hs{k) being indentically zero. Actually, if we assume a smooth behaviour as the quark 
masses go to zero, we can argue that the first class, corresponding to a purely scalar prop- 
agator {hi = /i2 = ^3 = 0), is empty. In fact, if one tries to construct such a solution by 
iteration, starting from the free propagator, the occurrence in the above limit of the strong 
infrared singularity due to the area law term (|1.4| ) is found to prevent the procedure to 
converge. 

As an elimination criterion of other classes one can use the unexpected occurrence of zero 
mass bound states in the corresponding BS equations and the consistency with Goldstone's 
theorem. Indeed in field theory the occurence of zero mass bound states appears as an 
exceptional circustance, expected only in connection with a spontaneously broken symmetry. 
In this way also the second {h-^ = 0) and the third classes {hi = /12 = 0) can be rejected as 
spurious and we are left with the fourth class alone, for which hQ, . . .h^ are all different from 
zero. Inside this fourth class various possible solutions can still be considered, which reduce 
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to the above classes when the quark masses mi and m2 vanish. By similar arguments also 
such solutions can be eliminated but one, that for which hi and /i2 keep different from zero, 
while /i3 vanishes for mi, m2 — > 0. 

The last solution corresponds to a breaking of the chiral symmetry in the vanishing quark 
mass limit and to the occurrence of a zero mass pseudoscalar solution of the BS equation, 
as one expects. Notice that what rules out the possibility of a solution of the first class, 
or reducing to the first class for mi,m2 —>■ 0, is the already mentioned infrared singularity. 
Therefore this is ultimately responsible for the chiral symmetry breaking. 

The plan of the remaining part of the paper is the following: In Sec. II, besides the first 
and the second order qq Green functions G^\xi,X2]yi,y2) and H^\xi,X2;yi,y2), already 
considered in we introduce the corresponding single quark propagators G^\x — y) and 
H^^[x — y). In Sec. Ill we discuss Eq. (|1.4P; we also introduce the functions H{xi, X2; yi, 2/2) 
and H{x—y), which are obtained by neglecting certain contributions related to the Schwinger 
strings, but become identical to H^\xi, X2] 2/1,1/2) and H^\x — y) in the limits X2 — > Xi and 
2/2 yi (notice that only such limits are relevants for bound states and chiral symmetry 
breaking). In Sec. IV we derive a SD equation for H{x — y) along lines similar to those 
followed to obtain the BS equation for H{xi,X2]yi,y2) in @1- In Sec. V we rewrite the BS 
equation in a form which makes it easier to compare to the SD equation. In Sec. VI we 
discuss the possible solutions of the SD equation and, in the limit of vanishing quark masses, 
the chiral symmetry breaking and the BS equation for a zero mass bound state. Finally in 
Sec, VII we summarize the results and make some additional remarks. 

Notice that no attempt of a numerical explicit resolution of the resulting SD equation is 
made in this paper. 



The quark- ant iquark and the single quark gauge invariant Green functions are defined 

as 



where -0'^ denotes the charge-conjugate fields, the tilde and Trc the transposition and the 
trace respectively over the color indices alone, U the path-ordered gauge string (Schwinger 
string) 



II. GREEN FUNCTIONS AND FEYNMANN-SCHWINGER 

REPRESENTATIONS 



G^\xi,X2, 2/1, 2/2) = -(0|T^^(a;2)^7(x2, xi)V'i(a;i)V^i(2/i)f/(2/i, 2/2)^2(2/2) |0) 
= ^Trc(?7(x2, 2;i)5'i(xi, 2/1; A)f/(2/i, 2/2)5*2(2/2, 3^2; -A)) 



(2.1) 



and 



G^\x - y) = {0\TU{y,x)ij{x)^Piy)\0) 
-iTrc(t/(2/,x)5(x,2/;A)) 



(2.2) 




(2.3) 
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(the integration being along the straight hne joining a to b, although otherwise specified, and 
P being the ordering operator on the color matrices), S, Si and 5*2 the quark propagators 
in the external gauge field A'^ and obviously 

\J[^\/ mA]Mj{A)e'SiA] ' y^-^) 

S[A] being the pure gauge field action and Mf{A) the determinant resulting from the ex- 
plicit integration on the fermionic fields (in practice however Mf{A) = 1 in the adopted 
approximation). Notice that we shall suppress indices specifying the quarks as a rule when 
dealing with single quark quantities. 

The propagator S is supposed to be defined by the equation 

(^7^D^ - m)S{x, y; A) = 6\x - y) (2.5) 

and the appropriate boundary conditions. 

As in [H we introduce the second order propagator 

{D,D^ + m^-]^g a^^F,,)A-{x, y; A) = -5\x - y) (2.6) 

and write 

S{x, y; A) = {i^D, + m)/V{x, y; A), (2.7) 

with {a^'^ = f [7^,7^^])- Then, after replacing (pTT]) in ( |2.2| ) and ( p.l|) , using an appropriate 
derivative it is possible to take the differential operator out of the angle brackets and write 

G^\xi,X2;yi,y2) = -{hidif, + mi){i-f2d2u + m2)H^\xi, X2;yi,y2) (2.8) 

and 

G^\x-y) = {i-f'^d^ + m)H^\x-y) (2.9) 

with 

H^\xi,X2;yi,y2) = -^Trc(^(x2, Xi)A^(xi, yi; A)[/(2/i, t/2)A^(a;2, 1/2; -i))- (2.10) 

and 

H^\x-y) = zTTc{U{y,x)A''{x,y;A)). (2.11) 

In the last equations 9^ stands for the ordinary derivative when acting on the propagators, 
while it refers to an appropriate distortion of the path when acting on a Schwinger string. 
In more explicit terms one can also write 

d,M{b,a)=d,M{b,a)- [\\\{bP-aP) \ r,U{b,a) (2.12) 

Jo dS'^P{a + X[b — a)j 
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(with SS^'^ = dz^Sz" — dz'^Sz'^ and the functional derivative being defined through an arbi- 
trary deformation, z ^ z + 6z, of the fine connecting a to 6 starting from the straight fine) 
and a similar expression for Squ. [] 

For the second order propagator we have the Feynman-Schwinger representation [ffl 



A''{x,y;A) 
with 



2 Jo 



rx rs ^ rs 

ds Vzexp[-i dT-{m^ + z'^)]S^F exp[ig drzf'AJz)], (2.13) 

Jy Jo 2 JO 



Sq = Texp 



A Jo 



- / dra 



SS^"'{z) 



(2.14) 



T and P being the ordering prescriptions along the path acting on the spin and on the 
color matrices respectively. Replacing ( ^.14| ) in ( |2.10D and ( ^.11[ ) we obtain 



and 



H^\xi,X2;yi,y2)= {-f 







dsi / ds2 







yi 



Vzi / VZ2 



y2 



I 

2 Jo 



dTi{ml + zl)- - / dT2{ml + zl)} 
I Jo 



exp I 

x5o^^5o^^^(TrPexp{^U dz^A^{z)}), 



(2.15) 



m\x - y)-- 



2 Jo 



ds / Vzexp { — - / dr^m^ + i^) 
2 JO 



x5o^(TrPexp{2(7i' dz''A^{z)}). 



(2.16) 



Here, the loop T^q occuring in the 4-points function is made by the quark world line 71, 
the antiquark world line 72 followed in the reverse direction, and the two straight lines X1X2 
and ^22/1 (Fig. as we already mentioned, on the contrary the loop Vq occurring in the 
2-points function is made simply by the quark trajectory 7 connecting y to x and by the 
straight line yx (Fig. 



^ Given a functional ^[7a6] of the curve 7a6 with ends a and 6, let us assume that the variation 
of $ consequent to an infinitesimal modification of the curve 7 — > 7 + ^7 can be expressed as the 
sum of various terms proportional respectively to 5a, to 5h and to the single elements 5SP^{x) 
of the surface swept by the curve. Then, the derivatives daP, d^p and 5/5SP^{x) are defined 
by the equation 5^ = da^daP^ + Sb^dbP^ + ^ SSP"" {x)5^ / 53^'' {x). For a Schwinger string we 
have 5U{b,a) = 6bPigAp{b)U{b,a) - 5aPU{b,a)igAp{a) + jl5SP"{z)F{-Fp„{z)U{b,a) and so 
daPU = -igUApia), B^pU = igAp{b)U and js^.U = V[-igFp„{z)U]. 
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III. WILSON LOOP WITH BACKWARD TRAJECTORIES 



Now we want to apply Eq. ( |1.2D to the evaluation of the Wilson loop integral occurring 
both in ( g7[5|) and in (|2:T6|) . 

At the lowest order the perturbative term can be written 

i{\n iy)pert = -\9^j dz^ / dz'^'D^^iz - z') (3.1) 

D ^u{z — z') being the free gauge propagator, which can be grafically represented by a waving 
line connecting the point z and z' on the loop. If we neglect the contribution coming from 
lines connecting a point on a trajectory to a point on a string or two point on the strings, 
we can write for Vqq 

4 /-Sl fS2 



4 /-si rs2 

ii\n.W)pevt = 7;9 dri / dT2D^^{zi - ^2)^1^2 - 

6 Jo Jo 

r dn r dr[D^,{z, - z[)z^z[' - ^g" C dr^ H dT^D,,{z2 - z'^)z^z',- + ... (3.2) 
O JO JO o Jo JO 



and for Tg 



z(lniy)p,rt = -^/ rdr r dT'D^,{z-z')z^z"' + ... (3.3) 
3 Jo JO 

Obviously the above approximation would not make sense in general, but we are even- 
tually interested in the two limits y ^ x and 1/2 ^ ^2, yi Xi and in which Eq.'s ( |3.3| ) and 
( p.2| ) become exact. 

Let us pass to consider the area term in ( |1.2D and refer to the loop Tgg first. 

In general one can writes 



S'min = / dt dX 



^ dt dt dX dx'^^ dt dx' \ ' ^ ^ ^ 

j^fj. _ ^iM(_x^i) being the equation of the minimal surface with contour Tqg. Since ( |3.4|) is 
invariant under reparametrization, a priori the parameter t can be everything. So we can 
assume that for fixed t one has 

«^(l,t) = ^r(n(t)), ti^(0,t) = ^2^(r2(t)). (3.5) 

However, if 71 and 72 never go backwards in time, we can also assume t to be the ordinary 
time, u^{s,t) = t. Then Ti(t) and T2{t) are specified by the equation 

z'Ari) = 4ir2) = t . (3.6) 

and the equal time straight line approximation consits in setting 

u%X, t)=t, M^(A, t) = Azf (ri(t)) + (1 - A)4(T2(t)) . (3.7) 

Replacing this in (|3.4| ) and performing a change of integration variables, one obtains ( |1.3| ). 
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As regards to Eq. ( |1.4j) , notice first that, as we already mentioned, this equation is 
identical to (|1.3|) , if the trajectories never go backwards in time. In order to discuss the 
situation we have when this is not the case, we assume that the extreme points coincide 
(xi = X2 and yi = 1/2) and the trajectories lie on a plane. Then let us refer to Fig. |^ 
and consider first the interaction term. For Z2 between the points R and S we have three 
different values of ti which solve ( |3.6|) and therefore three different intervals of ri which give 
contribution; these correspond to zi between A and B or between B and C or between C 
and D. In absolute value such contributions equal respectively the areas RABSR, RCBSR 
and RCDSR. However, when attention is paid to the sign factors appearing in front of the 
integral in A in ( |1.4| ), one realizes that the second contribution is negative, while the other 
two are positive. In this case the algebraic sum of the three areas equals simply the sum 
of CABC and RCDSR, and this amounts to the total area enclosed in the loop in the strip 
between the lines RA and SB. 

In a similar way, we can see that the contribution to the above strip coming from the 
self-energy term of particle 1 is also made by the sum of three terms (corresponding to z[ 
between A a.n B and zi betwen B and C or z[ between A and B and zi between C and D 
or z[ between B and C and zi between C and D), but in this case such sum is zero. 

An identical situation occurs for the strip between LE and xF. On the contrary, it is 
clear that no contribution comes from the interaction term to the area HFGH corresponding 
to a time larger than = X2, while this area is fully taken into account by the self-energy 
term. 

In conclusion, for Tqq on a plane and with the extremes of 71 and 72 coinciding, the Eq. 
( |1.4| ) is exact. For reason of symmetry the same must be true for a Tqq obtained from the 
preceding one by screwing the trajectories to elics along a line parallel to the time axis. In 
a general case ( |1.4| ) can be only an approximation, but we expect to be a good one, since 
the contribution coming from trajectories very far from the above classes should be strongly 
suppressed in the path integral (being far also from the minimun of the exponent). 

In a similar way in the case of we shall set 



•2 •/2 



Smm = -J^ dr dT'6{zo- z'q)\z- z'\e{zo)e{zQ) dX{: 

-(AfT4 + (l-A)iT^o)'}^ (3.8) 

and notice that even this equation becomes exact for y = x if Tg is of one of the types 
discussed above. 

As we shall see, for zero quark masses ( |1.4|) and (|3.8|) , taken together, are approximations 
consistent with chiral symmetry breaking and Goldstone theorem. 

Replacing (|3^ ) and (|L^ in ( p.l5| ) and (|0| ) and (|378D in ( p. 161) we obtain the following 
equations 

]^ poo poo pXl pX2 , i pSj 

H{xi,X2;yi,y2) = {-f dsi ds2 Vzi Dz2exp<^--V/ dTj{rin?- + z])\ 

Z Jo Jo Jyi Jy2 ^ J J ) 

2 

xSq^Sq" exp j dTj / dr'^Ei^Zj - z!^\ zj, z'j) - 

j=i 

-i dn dT2E{zi - Z2;zi,Z2)j (3.9) 
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and 



H{x — y)=- / ds I D^expj / (ir(m^ + i^)} 

2 Jo Jy 2 JO 

xS^ exp {i r r E{z - z'] i, z')} , (3.10) 
Jo Jo 



where we have set 

E{C-p,p') = E^,,,{C-p,p') + E,,^,{C.P,P') (3.11) 

with 

r ^pert =47r|«,D^,(C)pV^ 

I i^conf = 5(Co)|Ck(Po)e(p[,) /o' rfA{pgp[? - {Xp',p^ + (1 - A)poP'T]n^ ^ ^ 



Notice that in principle the quantities H{xi,X2; yi, 2/2) and H{x — y) as defined by ( p.9|) - 
( p.l2| ) may differ significantly from the original H^\xi, X2',yi, ^2) and H^\x — y) for arbitrary 
arguments, since, in writing such equations, we have neglected the contributions to the Wison 
loops coming from the strips betwen and X2, and ?/2, y^ and x^, for what concerns 
the area parts, and from propagators string-string or string-trajectory for the perturbative 
parts. However, as we mentioned, the two couples of quantities must concide in the limits 
X2 Xi, y2 yi, and y x. 

We can also define the quantities 

G{xi,X2;yi,y2) = -(^Tf'^iM + "ii)(^72^2v + ^-2)^^(2:1, X2; yi, 2/2) , (3.13) 

and 

G{x-y) =i{i-f>'d^ + m)H{x-y) , (3.14) 

where we have neglected the integrals along the string occurring in ( |2.12| ). Again 
G{xi,X2; yi, 2/2) coincides with G^\xi, X2; yi, 2/2) and G{x — y) with G^\x — y) in the consid- 
ered limit. 

Eq.'s (|3.9| ) and ( p.lO| ) are given in terms of purely configurational path integrals. For 
the need of the following sections path integral representations on the phase space are more 
convenient. These can be obtained from the preceding ones by performing a Legendre 
transformation on the exponents (see [Q for details). At the first order in E, or equivalently 
at the first order in cvg and ao? (a being a caractheristic lenght of the problem), we have 



00 



H{xi,X2]yi,y2) = {^Y I dsi I ds2 



/ VziVpi Vz2Vp2exp{iJ2 dTj[-ZjPj + -{p^j - rrij)]} 

Jyi Jy2 Jo I 

xS^'S^' exp {^Y.jj ^ ' dr'jEizj - z'j;pj,p'j) - 

-i dri dT2E{zi - Z2;pi,P2)} (3.15) 

JO JO 
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and 



]^ ;>oo rx rs ]^ 

H(x — y)=- / ds T)zT)pexp{i / dT[—zp-\ — (p^ — "W-^)]} 
2 Jo Jy Jo 2 

X exp {i dr'E^z ~ z';p,p')} , (3.16) 

JO JO 

which is justified by the fact that integration of the momenta in the gaussian approximation 
reproduces ^]9| and p.lO| up to the mentioned order. 

IV. THE SCHWINGER DYSON EQUATION 



From (|3.16|) a Schwinger Dyson equation can be derived by a technique strictly similar 



to that used in |^ to obtain a Bethe Salpeter equation. 

Using a trivial identity we can write with obvious notations 



H{x -y) = - ds VzVp exp [i / drK] 1 + 

2 JO Jy JO I 

+i dr I dT'E{z — z';p,p')exp[i df df' E{z — z;p,p')]> = 

JO JO JO JO J 

^ rOQ rx rs 

= Hq(x — y) -\ — / ds VzVp exp [i / drK] 
2 Jo Jy Jo 

rs rr rr rf 

X dr dT'Jab{z — z';p,p')a°'Sl,{exp[i df dTE{z — z;p,p)]} 

Jo Jo Jt' Jt' 

xcr^iSg {exp[z / df df' E{z — z';p,p')] 
Jo Jo 

X exp [i r df r df'E{z - z';p,p')]} , (4.1) 
Jt' Jo 



where 



and 



K = -pz + p^ -m^ (4.2) 



^ /* OO px ps 

Ho{x-y) = - ds VzVpexp[i drK], (4.3) 

2 Jo Jy Jo 

We have also set a,b = 0, fiiy, with = 1, and have denoted by Jab coefficients which come 
from E and its commutation with Sq and are given by (see Eq. (A8)) 

Jo;o = E{z - z';p,p') = 

= An^a^D^fsiz - z')py^ + 

rl 1 

+a5{zo - 4) \z- ^\ e(po)e(Po) / dX{plp'^ - [Ap'oPx + (1 - \)voV'vfY 

Jo 



4 



Jp.v,o = -7T^as{S"du - 6^d^)Daf3{z - z')p' 



3 
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-a6{zo - ZQ)e{po) —^^—-== 



Jo;pa = vr-asP"(5;5. - S:d,)D^p{z - z') + 

+a5{zQ - Zq)po — i^— ^==== — ^e(po) 

J;..;p. = -7^«,((5^a, - - - /) (4.4) 



Up to the first order in E we can replace the last exponential in ([4.1| ) by 1. Then we 
obtain 

H{x -y) = Ho{x -y) + i J d^id^rid^i' Hq{x - 

xU{i,i';v.V>''H{r^-r^')a'H{i' -y) , (4.5) 

where 

d'p d'p\-^p(^-r,)J^^^i±V_e+^.p^p>^^-^p'(,'-e) (4.6) 



(27r)4 (27r)4 
and we have used the equation 

/; VzVp exp [t drK] ... = J d^id'r^ J (0e-^(«-'') 

X J^VzVpexp[iJ^dT'K]J^VzVpexp[iJf;dT'K]... , ^ ' ^ 

which follows immediately from the discrete form of the functional integral. 
Furthermore, by performing a Fourier trasform of Eq. ( |4.5|) we obtain 



H{k) = Hoik) + 

+^ / Mk) L{k - I; ^) a'H{k)- (4-8) 

Finally, by comparison with 

H{k) = Hoik) + iHo{k)t{k)H{k) , (4.9) 
]8[) gives the Schwinger Dyson equation 

f (A:) = j^M^- ^' • (4-10) 
In an explicit form the kernels lab are given by 
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+ 4a I d'Ce''^<\C\e{po)eip',) d\{plp'^ - [\p',p^ + (1 - A)pop't]'}^ 



U,pa{Q-^P,p') = -47rz-a3P"(5?g. - 5fg,)D«^(Q) + 
where in the second and in the third equation Co = has to be understood. 



(4.11) 



Notice that, were |4.10| covariant, i.e. were the quantities lah tensors, we could write 

iH-\k) = a(A;2)+ /:b{k'^) . (4.12) 

However, due to the equal time straight line approximation, that is not the case; we must 
assume ([4.11 ) to hold in the meson center of mass frame and set 



r=0 



with 



Wo = 1 , = 7° , i02 = -j-k, ^3 = 7°7 ■ A; , 



(4.13) 



(4.14) 



k = -r^k and /io(A:), . . . h-^{k) functions of ko and \k\ separately. Obviously ( |4.13|) reduces to 

\k\ 



(|jg) for ho = a{k'^), hi = ko 6(fc2), h2 = \k\ b{k'^), h^ = 0. 
From (pD and (|11) it follows 

3 

f(A;) = iH-\k) - iHQ\k) = ho{k) - - m^) + Y.Ur{k)hr{k) 



r=l 



and 



H{k) = i- 



ho{k)-Y.l=iOJr{k)K{k) 



^ hl{k) - hl{k) + hl{k) - hl{k) ' 
Replacing such equations in ( [4.10|) and taking advantage of 



we finally obtain 



(4.15) 



(4.16) 



(4.17) 
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hr{k) = Sro{k'^ — nn?) 



3 
s=0 



dH 



Rrs{k,l) hs{l) 



{2nYhl{l)-hl{l) + hl{l)-hl{l)' 



with 



Rrs{k, I) = T\lai>{k -l;^, ^)TT[u;nkKus{k)a'] 



(4.18) 



(4.19) 



where the sign - apphes to the s = case, the sign + to all the other cases. Notice that 
from well known properties of Dirac's matrices it follows immediately Rqi = Riq = Rq2 = 
R20 = Ri3 = R31 = R23 = R32 = 0, while from 



^Mr.;0(Q;g,g) 



{Q;q,q) and {Q;Q,q) {Q;Q,q) 



(4.20) 



it follows also i?03 = -R30 = 0. Actually only Rqq, Ru, R12, R21, R22, R33 are different from 
zero. Expressions for these quantities in terms of the lab are given in Appendix B. Here, as 
an example, we want only to report the explicit expression 



Roo{k,l) = -47r|as[4 



(p-/)4 ^ 



(4.21) 



-a/ Ae-*('-')-^|CII(^o + /o)lV(^o + - {kT + h? 
The confining term in this equation contains an infrared singularity essentially of the type 

^2 1 



de^ {k - ly + e2 



(4.22) 



with e — > 0. Similar singularities appear in the other diagonal kernels Rrr- 

Expressions like ( |4.22| ) correspond to well defined distributions and are harmless if they 
appear in integrals multiplied by regular functions. However, they can give troubles and 
generate new singularities when occurring in connection with discontinuous functions. As 
it appears from (^4.21|) and ( [4.11|) , R^s have also a bad ultraviolet behaviour that should be 
handled by renormalization. Renormalization would affect the perturbative parts of R^s-, 
but its explicit consideration is not important for our present purpose. So we shall simply 
suppose to regulate ( |4.18|) by a cut-off A. 

Finally going back to the first order propagator ( p.l4| ) we can also set 



iG-\k) = Y,ujr{k)g,{k) 



(4.23) 



and then we have 



r=0 



' hQ = mgo + kogi + \k\g2 
hi = mgi + kopQ - \k\g-:i 
h2 = mg2 -\k\go + k'^gs 
h = mg^ + \k\gi + kog2 



(4.24) 



From ( [4. 231 ) it is apparent that, for zero quark masses, chiral symmetry breaking corre- 
sponds to go{k) 7^ and/or (73 (/c) 7^ and consequently from ( [4.24| ) to hi{k) and h2{k) not 
simultaneously vanishing. 
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V. THE BETHE-SALPETER EQUATION 

The Bethe-Salpeter equation for the quantity H{xi,X2;y 1,112) is derived from ( p.9|) by 
applying the same identity used in ( |4.1| ) to the interaction, rather than to the selfenergy 
term, as done in the case of the Schwinger- Dyson equation. In terms of the quantities defined 
by Eq.'s (|46|) and (U) it reads 

H{xi, X2] yi, 2/2) = Hi{xi - yi) H2{x2 - 2/2) - 

-i j d'^^id'^C2d\d%Hi{xi - ^1) H2ix2 - 6) 

x4fe(6,6;^i,??2Xo-2i:/'(?7i,?72;2/i,?/2) , (5.1) 

where Hi and H2 denote the quark and the antiquark propagators as defined by ( |3.1CI| ). 
The corresponding homogenuous equation in the momentum representation is 

L(k -n.\p + \p - :^)<aS*p(«) , (5.2) 

the center of mass frame being understood, i.e. P = (m^, 0). 

The wave function $p(/c) can be also reinterpreted as a matrix (in which the column 
index refers to the quark and the row one to the antiquark respectively) and the pedices 1 
and 2 in the spin operators suppressed. Then setting 

<l>p(fc) = i/i(ip + A;) T'j^ik) H2{ -^P + k)C (5.3) 

and using C H'^ {k)C^^ = H'^{—k) we can rewrite ( p.2| ) as 

(k) = -i / - — -lahik -u,-P+ , -P 

' J (27r)4 ""^ ' 2 2 ' 2 2 ' 

cr- H^{lp + u) r^^(«) H^{-]^P + u) (5.4) 

with 

a° = aV a'^" = Cia^'^'fC-^ = -a"" , (5.5) 

C being the charge conjugation matrix. This is the form of the BS-equation more convenient 
for our present purposes. 

More specifically we need to consider the case P = 0, corresponding to a zero mass 
bound state. If we take into account the property (cf. ( |4.11| )) 

iab{Q;p,p')(T'' = -iabiQ;p,-p'y, (5.6) 

and assume mi = m2 = m, we have 

r^ik) = ^ j - ^' ^-^)a^H{un{u)H{u)a' . (5.7) 

Notice that a zero bound state exists, if (|5.7|) has a non trivial solution. 
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VI. CHIRAL SYMMETRY BREAKING AN GOLDSTONE THEOREM 



To discuss Eq.'s ( |4.1CI| ) or( ^.18| ) and ( |5.7D it is convenient to rewrite them in euclidean 
form. This can be obtained by making the substitutions po —>■ ip^ in all momentum variables 
and by setting 70 = 274. To solve the ambiguities related to the occurrence sign funtions and 
the square roots in ([4 .111) one should go back to the eucleadean counterparts of the original 
equations and particularly of ( p. 15 ) and (|3.16|) . Alternatively one can extract from the 
square roots factors as \po\ and then set e(po)bo| = Po before performing the substitution. 
It also convenient to redefine hi as ihi, /13 as ih^ and also R12 as iRi2 and R21 as —iR2i- 
The resulting expression for Rqq is, e.g., 



Roo{k,l) =-47r|a,[4 



(p-iY 



+ iJ + 



|C||(A:4 + /4)IV(^4 + /4)2 + (A;t + /t) 



(6.1) 



while expressions for the confinement parts of the other Rrs are reported in appendix C. 



Then, making explicit the terms that are different from zero, (^4.18|) can be rewritten 



( U - _p _ ^2 I r dH Rooik,l) hojl) 



h (h\ — r d.'^l Rii{k,l)hi{l)+Ri2{k,l)h2(l) 

"-lyf^) - J (27r)4 

h (h\ — f d-^l R2i{k,l)hi{l)+R22{k,l)h2{l) 

r^yti) - J (2^)4 h'^^(i)+h-i{i)+h-i{i)+hi{i) 



h (t.\ _ f dH R33(k,l) hsjl) 

Similarly for the euclidean counterpart of 



(6.2) 



3(0 

we have 



d'^u f k + u k + u, ho{u) + ^Jii{u) + ^ ■ uh2{u) + ■ uh'i{u) 

Ub{k - u, , )a — 



{2tt 



Ml ^ ho{u) + 74/^1 (m) + 7-M/i2(n) +747-^^3(m) 5 
° ^""^ hliu) + hl{u) + hl{u) + hl{u) 



(6.3) 



Notice that a priori Eq. (|6.2| ) could be expected to admit various types of solutions 
corresponding to different sets of identically vanishing 



Al) 
A2) 
A3) 
A4) 



/ll = /I2 = /l3 = , 

hi, /i2 7^ 0, /13 = , 
/ii = ^2 = 0, /l3 ^ , 
hi, h2, h-s^O. 



(6.4) 



To gain more insight on the nature of such solutions let us assume that we can solve 
2]) by iteration starting by an initial ansatz of the type 



(0) 



k. 



m 



hf^ = ^lu-) 



1,2,3) (6.5) 



/ii,/i2,/U3 and fL being constants with the dimensions of a mass and fr arbitrary functions 
chosen in such a way that at least the first iteration is meaningfull. E. g. after the first 
iteration we have for ho 
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"■n (k) = —k — m — - — — t^t 777^ jr-. . (6.6) 

° ^ ^ J i27T)^(P + m''y + [hf\l)Y + [hPil)f + [h'i\l)]^ ' 



Solutions of the types A1-A4) can be obtained from ( |6.5|) simply by taking /ii = /X2 = 
/i3 = 0, //i, /i2 7^ 0, /i3 = etc. respectively. Notice that fii, . . . p,, when different from zero 
must be expressed in terms of the masses existing in the theory, m, y/a and the cut-off A 
or the renormalization scale. 

Let us consider separately the various cases. 

Case Al). In this case we expect no actual solution smooth in m = 0. In fact, setting 
m = and hi = h2 = = Eq. (|6.6|) becomes 



h^^\k) = -kl -P- [ fl^^oikj) . (6.7) 

Having in mind Eq. ( ^4.21| ) we realize that integration over I4 in such an equation produces 
a singularity in the tridimensional integral which conspires with ( [4.22|) for k ^ bearing 



a singularity in h'^\k). Similar circustances occurre at various steps in the iteration 
rendering the entire process instable and apparently preventing convergence. 
Case A2). One can immmediately check that, if vanishes, setting 

T^{k) = [jMk)+l-kh{k)]j,, (6.8) 

( |6.3D is made identical to the subsystem formed by the second and the third Eq. ( |6.2D and 
so non trivially verified. Therefore, a zero mass pseudoscalar meson exists, independently 
of the quark mass value and of the occurrence of a specific symmetry. This class has to be 
rejected. 

Case A3). This time a non trivial solution of (|6.3| ) can be obtained by setting 

T^'ik)=-fMk)- (6.9) 

This would correspond to a zero mass scalar meson. Also this class must be considered 
spurious. 

We are left with case A4 in which all the quantities hi, /i2 and h^ are different from zero. 
Inside this class we can still distinguish various, possibly single, solutions corresponding to 
different behaviours of the hs in the limit m — * 0. Schematically we can think of generating 
such different solutions by setting in (|6.5| ) 

/is = m, 

v^' ^=^=;i^' (6.10) 

Notice that Bl, B2, B3 correspond to solutions that are reduced respectively to the types 
Al, A2, A3 for m — *• 0; in fact the factors m in front of the various fr{k) are reproduced at 
every step of the iteration. Notice also that for m — ^ the theory acquires chiral symmetry. 
As we mentioned the symmetry is broken if in the same limit hi and ^2 do not simultaneously 



Bl) 


^1 


= 1^2 


B2) 


^1 


= /W2 


B3) 


/Wi 


= fJ'2 


B4) 


/Wi 


= fJ'2 
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vanish. So we have a broken chiral symmetry in cases B2 and B4, actual symmetry in the 
other two cases. 

Let us consider again the various cases. 

Ansatzes Bl and B3 reproduce similar situations as Al and A3 for m — and can 
be excluded for the same reasons. Ansatz B4 contradicts Goldstone's theorem; in fact for 
m ^ 0, /i3 keeps different from zero in this case and the expression (|6.8|) would not be a 
solution of ( |6.3D . On the contrary ansatz B2 seems to satisfy all requirements. For m — *■ 
( |6.8| ) becomes a non trivial solution of ( |6.3| ) and and a zero mass pseudo scalar bound state 
occurs. 

In conclusion we are left with solution B2 alone. This seems to be the only one which 
is physically sensible and mathematically consistent. For light u and d quarks it would 
correspond to a breaking of the approximate chiral symmetry and would correctly provide 
a TT meson with a small mass. 



VII. CONCLUSIONS 



Using the formalism developped in 0] , the minimal area law and the straight line approx- 
imation in a modified form, we have obtained the Schwinger-Dyson equation ( [4.1(J| ) for the 
quark propagator occurring in the Bethe-Salpeter equation ( |5.1| ) or ( |5.4| ). We have explicitly 



rewritten the SD equation as a system of non linear integral equations involving four scalar 
functions. However, we found that such system has not a unique solution and had to resort 
to additional criteria to select the correct one. 

After a first selection, assuming an iterative resolution procedure to converge, we have 
been brought to consider various possible solutions corresponding to different starting 
ansatzes. One of these (solution B2) produces chiral symmetry breaking in the quark mass 
vanishing limit and a zero mass pseudoscalar meson. This corresponds to what we believe 
to be the real situation. All the other solutions are mathematically inconsitent, or in contra- 
diction with general theorems, or, simply, have unexpected features which we do not believe 
should be present in real QCD. If we accept B2 to be the correct solution, then, we achieve 
a very consistent formalism for the treatment of the quark-antiquark bound states in QCD, 
both for light and for heavy quarks. 

We stress once more that the occurrence in the kernels of the infrared singularity ( |4.22| ) 
is what seems to prevent solutions Al or Bl form existing. Therefore a strict connection 
appears to exist between confinement as expressed by the area law and chiral symmetry 
breaking. 

Notice, also, that the possibility of a solution with the properties we have found depends 
strictly on the assumption of ( |1.4| ) in place of (|1.3|) . Had we used straight line approximation 
in the original form, only the perturbative parts of ( 4.11 ) would have occurred in ( |4.10D . 



Then the kernels in ( |5.7| ) would have been different from the kernels in ( [4.11| ) and 



would not have been a solution of (|5.7| ). The correct solution of (|6.2| ) would have been of 
the type Al and chiral symmetry breaking and the zero mass pseudo scalar meson would 
not have occurred. 

Let us make a final comment. 

As we mentioned the confinement parts of the kernels lab, as given by ( |4.11| ), are not 
tensors. As we mentioned, this circustance is a consequence of the equal time straight line 
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approximation, which is obviously frame dependent. We have assumed the priviledged frame 
to be the center of mass frame of the meson. But then lab can be rewritten as tensors, if we 
introduce exphcitly the total four-momentum of the meson P. For this aim it is sufficient 
to set 

Q\\ = Q'l = Q'~ Q\\^ (7-1) 

give analogous definitions for p and p' and make the identifications 

Ql^Ql Q'^-Qi etc. (7.2) 

and similar ones. 

In the above perspective the quark propagator H{k) should be rewritten as H{P, k) and 
the Eq. (|4l3| ) as 

iH-\P, k) = j2 Lo',{P, k)K{P, k) , (7.3) 



with 



r=0 



J^ = l. uj\=^P, J^= (7.4) 



and ( |6.2| ) correspondingly modified. 

Since the quark propagator has been originally defined in an independent way by ( p.llD 



at first sight it can be surprising that now it is made dependent on the meson momentum. 
Notice, however, that it aquires again P dependence as a consequence of the straight line 
approximation. After that it becomes the specific propagator to be used in that BS equation. 

Notice also that the need to choose solution B2 (with /i3 7^ for m 7^ 0) is a consequence 
of the above dependence. Had lah^ been tensors in k and / alone (as in the abelian case), we 
should have had /is = (see Eq. (4.12)) and Al would have been the only possible sensible 



solution. Therefore a P dependence of the quark propagator seems to be essential for a 
consistent chiral symmetry breaking and not an artifact of the minimal area law model or 
of the straight line approximation. 
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APPENDIX A: FUNCTIONAL DERIVATIVES 

Taking into account appendix C of Q], and specifically Eq. (CIO), we can write for the 
perturbative part 
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6S''''{z{T))Ja 



df I df'Epcrtiz- z',z,z') 



4vr^as{x(a,6)(r) £ df {5^^d, - 5:d,)D^ 



ai3{Z — z')z^ 



-Xia',v){r) / df{5ld, - 5^,d,)D^p{-z - (Al) 



where X(a,b){.T) and X(a',b'){T) are the characteristic functions of the intervals (a, h) and (a', 6') 
{{a,b) and {a',b') contained in (0,s)) and the derivatives are intended to act on the entire 
argument of -Da/j. Due to the arbitrarity of the above intervals, we can also write formally 

<5Sf"'(2(r)) ^pert{z — Z ,Z,z) = 

= AnIaMr - - 6'^d^)D^p{z - z')^'^- (A2) 

-Sir - f')(5^a. - S^MDa^i^ - z'V] . 

For the confining part for consistency we must define 

5 



5S^'^{z{T)) 



df df'Econi{z - z,^,z) (A3) 

J a J a' 



as the the straight hne approximation of j§jn;S^^^,, the quantity S^J^^, denoting the area of the 
strip of the minimal area enclosed by Tq delimited by the two geodetics connecting z{a) to 
z{a') and z{h) to z{h') respectively. Then, having again in mind Eq. (C9) of 0, we have 

(5SMf (^(r)) ia Sa' df'Econiiz — z',Z,z) = 

■ {z'^-z''^)z''-{z''-z'")z'^ 



(^X{a,b)ir)X(a',b')ir')€{zo)€{z'Q) 

{zt^-z'i^)z'" -{z'' -z"')z'i^ 



z'2(z-z')2+{i{z-z'))^]i 



[-i'2(2-z')2 + (i'(z-2'))2]i 



2()(-r')=2o(T) 



(A4) 



a\x(a,b){'^)Ia' dr'S{zo - ^o)|4|e(^o)e(4) !Z/2l\ .^f-. 52a 

I, lz^{z-z )^ + {z{z-z 

-X{a',b')iT) !adf5{zo - zo)\^o\e{^o)e{zo) 



b^^Xf^. _ ^ (z^-znz-^-jz-^-znz^^ 

[i2(2-5')2 + (i(f-2))2]2 



or 



Econii^Z Z , Z, Z 



= a6{T - f)6{zo - ZQ)e{zo)z'Q 



-a6{T - r')6{zo - 4)e(io)i;o 



^/ {z^^■~z''')2''-{z•'-z'•')2'' 
2 -J 

Z-—Z I \ / -2^0 — 



(A5) 



1^ ^1, /-■ /2 i'2 

\z-z ly 20 — 



Notice also 



= -8^|aJ5(r - f)b{T' - f'){bld, - bZd^d, - 6%)D^p{z - z') ^ ' 
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and 



Then 



SSp''{z{T'))5St'''{z{T)) 



^coniiz Z , Z 



lS^,E{z - z\ i, z'){S',)-' = E~ \a^^ f 

4 JO 

+— a'^v" / de / de' 

32 io io 



5SP<'{C{0')) 



de 

6 



ssp'^iCie)) 



E 



E 



+ 



E{z - /, i, z') - i(T^'^|47r^as(5^<9, - 5:d^)D^p{z - z')z'l' + 
+a5{zo - 4)e(^o)4^ / ^ ^ 

I ^ ^ I V ^0 ■^T 

-47r^as - d^MD^piz - z')z'' - 

, iz^ - z'f'^z"' - iz'' - z"')z''' 
-(t6[zq - Zf^)e[Zf^)zQ 



\Z — Z \ \l Zq — z. 



U2 1^2 
T 



4 
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[z — z 



(A7) 



(A8) 



Finally, after the Legendre transformation, z^ and zl^ are simply replaced in (|A8|) by p^^ 
and p''^. 



APPENDIX B: EXPLICIT EXPRESSION OF THE SD KERNELS 



Evaluating explicitly the trace in Eq. ( [4.19| ), one obtains 



Roo{k,l 
Ru{k,l 
Ri2{k,l 
R22{k,l 

[ R33ik,l 



-/o;o - 2gPPg-'^I^,,p^ 

/o;0 - ?>g'"'hu-fia + 2gt'Pg''''I„ 

-R2i{l, k) = Aik,J - Sg'^^iou-J' 

io-ik ■ i+8g''^ii,.jM' + 2gPPg^"i^,,p,k ■ / + 2z(/o;i,- - hy^^UV 
Ufi - 8ig'^ + 2FP)4;0, + Ag'^iiu.jM + 2gPPg''^ 1^,,,^ 



where the arguments are intended as in ( [4.18| ). 



(Bl) 



APPENDIX C: EUCLIDEAN FORM OF THE CONFINING KERNELS 



The euclidean form of Rrs{k, I) is obtained from Eq. ( |Bl| ) and ( [4.1 1| ) by setting p = p' = q 
and then making the substituitions Qo iQ^, Qq iq^, R12 iRu, R21 ^ —iR2i, 
according to the rules given in Sec. IV. 

We report the explicit expression for the confinement parts, the only non trivial: 
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R^,rik, l) = -aj d'Ce-'^<\q,m ^/dT^ + 4z ^^''i^^^^^'l^''^ ] 
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FIGURES 



FIG. 1. Wilson loop Tgg for the four points Green function formed by the quark and the 
antiquark trajectories (full lines) closed by two Swinger strings (broken lines). 

FIG. 2. Wilson loop Tq for the quark propagator formed by the quark trajectory closed by a 
Scwinger string. 

FIG. 3. Loop Tqq for conciding end points with the quark trajectory going also bakwards in 
time. 
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